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We present an extended study of finite width zigzag graphene ribbons (ZGRs) based on a tight- 
binding model with hard-wall boundary conditions. We provide an exact analytic solution that 
clarifies the origin of the predicted width-dependence on the conductance through junctions of 
ribbons with different widths. An analysis of the obtained solutions suggests a new description of 
ZGRs in terms of coupled chains. We pursue these ideas further by introducing a mapping between 
the ZGR model and the Hamiltonian for N coupled quantum chains as described in terms of 2N 
Majorana fermions. The proposed mapping preserves the dependence of ribbon properties on its 
width thus rendering metallic ribbons for N odd and zero-gap semiconductor ribbons for N even. 
Furthermore, it reveals a close connection between the low-energy properties of the ZGR model and 
a continuous family of square lattice model Hamiltonians with similar width-dependent properties 
that includes the 7r-flux and the trivial square lattice models. As a further extension, we show that 
this new description makes it possible to identify various aspects of the physics of graphene ribbons 
with those predicted by models of quantum spin chains (QSCs). 

PACS numbers: 75.10 Pq, , 3.20 At, 73.63.-b,71.30.+h ,73.43.f 



Graphene has attracted a lot of interest since its suc- 
cessful synthesis in 2004i. The crystalline structure of 
the material is given by a triangular Bravais lattice with 
a two-atom basis which provides a spinor character to 
wave-functions. Calculations of the band structure re- 
veal that conduction and valence bands coincide at six 
corners of the Brillouin zone (the so-called Dirac points). 
The linear spectrum around these points makes possible 
to describe the low energy physics by Dirac-type equa- 
tions. Therefore graphene shows many unique properties 
not observed in other materials. Among them we should 
cite the unusual quantum Hall effect and the high con- 
ductance of charge carriers through potential barriers, a 
property which is attributed to the Klein paradox^. 

From the point of view of potential applications, a 
high mobility and stiffness makes graphene an excel- 
lent candidate for future electronic devices. Interestingly 
enough, confined graphene samples (such as ribbons and 
quantum dots) reveal even more peculiar properties that 
strongly depend on edge terminations^. Thus, a cut- 
ting through armchair edges, renders ribbons that can 
be metallic or semiconducting, depending on the rib- 
bon width, while ZGRs have localized states near the 
edges with an almost flat band between any two in- 
equivalent Dirac points 3 -. Furthermore, numerical^ and 
theoretical^ work predicted a conductance of ZGRs 
strongly dependent on the ribbon's width. These works, 
focused on transport properties through junctions, con- 
cluded that junctions between ribbons with even and odd 
number of zigzag lines showed a valley- valve effect. 

In this work we show that this even-odd behavior has 
a simple explanation based on the presence or absence 
of zero modes in the spectrum of a finite width rib- 
bon. This realization leads naturally to a description of 
ZGRs in terms of coupled chain models that can be ex- 



tended (in some cases) to models of quantum spin chains 
(QSCs). This dual description suggests that ZGRs can 
be viewed as a material realization of many different the- 
oretical models of QSCs. At the same time, the well- 
known physics of various QSCs models can shed light on 
various properties of graphene ribbons. 



I. MODEL FOR ZIGZAG GRAPHENE RIBBONS 

We first review some results obtained with a tight- 
binding model for ZGRs with hard- wall boundary condi- 
tions. In the absence of magnetic fields and spin-orbit in- 
teractions, the Hamiltonian is spin-independent and thus 
it is enough to consider spin up (or down) electrons only. 
The Hamiltonian and the spinor wave-vector in momen- 
tum space are given by: 

\<p ) W± °^ B = ±e~ M / 2 J 

(!) 

with (fi(k x ,k y ) = tj_e tkyb + 2t cos an H <p(k x ,k y ) = 
<p{k X i ~k y )- Here k x , k y are measured from the center of 
the Brillouin zone (the T point) and b = av3/2. The 
angle a is defined by (p(k x , k y ) = y*(p(pe la . By an appro- 
priate gauge transformation, all atoms along each zigzag 
line can be labeled by the same coordinate y. This is 
equivalent to work with the deformed lattice shown in 
Fig[T] Notice that in graphene the inter- and intra- chain 
couplings arc equal = t). However for reasons to be 
clear later, we introduce these couplings as independent 
variables. 

Standard hard- wall boundary conditions for ZGRs are 
given by ua(v = —W/2) = and Ub{v = W/2) = along 
the lower and upper edge respectively- ^ 13 ' 14 . It is useful 
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FIG. 1: The boundary conditions for ZGR implies that ua = 
on the lower edge and ub = on the upper edge (green 
dotted line). All atoms alone each zigzag line are labeled by 
the same ^-coordinate. The right panel shows a generalized 
lattice model with interchain hoping term t± and extra hoping 
term r)t±. In graphene 77 = and t± = t. 



to introduce the variable N that counts the number of 
zigzag lines or chains inside the ribbon as N = W/b — 1. 
Note that due to translational symmetry along the re- 
direction, k x is a good quantum number, and for a given 
k x the wave function is given by the linear combination 
of two degenerate states at ±h y 



sm(k y y + a/2 — mr/2) 
sm(k y y — a/2 — rwr/2) 



in which k y satisfies 



k y W 



mr. 



(2) 



(3) 



where n is an integer. The condition imposed by Eq. © 
can also be written as tp/tp = c 2lkyW . This also implies 
that 



Ak v W 



sin k y b/ sin k y W. 



(4) 



One peculiar feature of ZGR is the presence of bands 
with complex wave number k y (and correspondingly com- 
plex a) between two in-equivalent Dirac points K = 
(27r/3, 0) and K' = (4-7r/3, 0). Previous works in the limit 
of semi-infinite ribbons have shown that these bands are 
flat and correspond to highly localized states at the edge 
of the ZGR~ & 13 i 14 . These results can also be easily ob- 
tained from the definition of ip and Eq. Q which supports 
a solution with e lkyb ps — 2(i/£j_) cos(k x a/2) with zero en- 
ergy. 

To get further insight into the properties of finite width 
ribbons we derive below the expressions for wavefunc- 
tions and energies of these quasi-degenerate bands. 

First we consider the case where k x > it that renders 
cos(k x a/2) < 0. For a zero energy solution to exist (ip = 
0), it is necessary that e zky > 0. This implies that k y 
should be imaginary and we can write k y = iq where q is 
a positive real number. By returning to Eq. ^ we obtain 
e M > which implies that a is imaginary. Therefore the 
condition set in Eq. (J3J) can be satisfied by setting n = 
and the wave- functions for the edge states are given by 



X±(k x a > 7r)) 



sinhq(y + W/2) , 
1 ±sinhg(y- W/2) ' ' [ ' 



Now consider the case for k x < tt for which 
cos(k x a/2) > 0. The zero energy condition now requires 



that e lkyb w —2cos{k x a/2)t/t^ < and thus k y must 
be a complex number: k y = tt + iq where q is a posi- 
tive real number. In this case Eq. §5§ can be satisfied by 
n = N = W/b — 1 which results in 

v+(k a < tt) - C ( ^[{q + ^){y + W/2) \ 

The energy dispersion for all bands is given by E = 
±t± sm(k y b) / sin(k y W) . In particular, the edge state dis- 
persion near k x a = n can be approximated as: 



E « ±t±(tk/tj 



(7) 



where k it the dimcnsionless momentum measured from 
the center of the band k x a = tt (located between two 
inequivalent Dirac points). 

Eq. ([7]) reveals an important feature of ZGR: ribbons 
with an odd number of chains possess different transport 
properties than those with an even number of chains^. 
This feature is reflected in the phase factor of Ntt that 
appears in Eq. ((SJ) which changes the relative sign of 
the spinor components in the wavefunction. A clear 
picture of this feature emerges when considering wave- 
functions at three different y— points (one at the center 
of the ribbon y = and two points on opposite edges 
y = ±W/2) and finding their asymptotic form when ap- 
proaching k x a = tt. 

The argument proceeds as follows: first, evaluate the 
wavefunctions at k x a > tt and obtain their limiting ex- 
pression as k x a — ► tt. 

At y = as q — > 00 the spinor is: 



X±(k x a > tt) ps C 



,qW/2 



qW/2 



CeqW /2 



1 



(8) 



C ~ e~ qW and the wave function inside the ribbon goes 
to zero at k x a = tt. Notice the relative sign of ua/ub- 

At the left edge y = —W/2 as q — > 00 the expression 
reduces to: 



c 







To 



qW 



Ce qW 







(9) 



and finally, at the right edge y = W/2, as q — > 00 we 
obtain 



xl 



c 



e qW 





Ce qW 



(10) 



Note that for each band (conduction or valence) the 
signs of the spinor components are consistent when we 
start from the left (Eq. ©) go through the middle 
(Eq. ©) and then to the right (Eq. (fTD])). 

Now let us approach tt from the left, i.e., with k x a < tt: 
At y = as q — > 00 the expression for the spinor is: 



X±(k x a < tt) ps C 



e qW/2 e iKW/2 
± e qW/2 e -nrW/2 



n qW/2 e ' 



(11) 

(12) 
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At the left edge y = —W/2 and as q — > oo 



the Hamiltonian becomes Hm = / dx[Ho + Ti'\ with 



xi 



c 





±ffl w 



Ce qW 





±1 



(13) 



2/V 



Similarly, at the right edge y — W/2 and as q — > oo 



c 







Ce 9W/ 



o 



(14) 



For W odd (AT even) any two conduction bands x+ , X+ 
have the same asymptotic wave function (notice the ra- 
tio of ua/ub = — 1). The same results holds for any two 
valence bands wavcfunctions with a ratio ua/ub = 1- 
This means that as k x moves continuously through the 
band-center, the conduction band does not cross the zero- 
energy (Fermi energy) point (the same holds for the va- 
lence band). In contrast, for ribbons of W even (N 
odd), two bands x+i X— with different energies are the 
ones that have the same asymptotic wavefunction (no- 
tice the ratio of ua/ub = — 1). Similarly, the two other 
bands \+ have their wavcfunctions with the ratio 
ua/ub = 1 as the center of the band is approached. This 
means that if we change k x from left to right we go from 
conducting band to valence band or from valence band 
to conducting band, i.e., the bands cross. 

From Eqs. (J5J) and ©, we see that for odd N the 
overlap between two spinors in the conduction band 
at k x a = 7r ± k is equal to zero: J dy < x+(k x a = 
ir+k)\x+(k x a = Tt—k) >= 0. On the contrary, for even N 
we have J dy < x+(k x a = tt + k)\x~ (k x a = ir — k) > = 0. 

The analysis of finite width ribbons presented here 
clearly shows that the origin of the even-odd dependence 
found in previous works in the limit of semi-infinite rib- 
bons, is given by the nature of the band of the zero energy 
mode. 



II. GRAPHENE RIBBONS AS COUPLED 
QUANTUM CHAINS 

The even-odd dependence formalized in the previous 
section suggests that some physical properties of ZGR 
should be captured by highly anisotropic models where 
a ribbon of finite width is viewed as a set of coupled one- 
dimensional chains. As shown in RefJ^, the Hamiltonian 
of N decoupled (fj_ = 0) chains is Hn = Y^n=i ^ n wnere 
H n = J dx\tcA n ]{x)(cB n {x) + CB n (x — a)) + h.c] is the 
Hamiltonian of the n th chain. The energy of each chain 
is given by E = 2tcosk x a/2. Near the center of the 
band k x a = n, we can write the lattice creation and 
annihilation operators by: 

c An (x) =i(-l) x / a (R n (x)+L n (x))/V2 
c Bn (x) = {-l) x ' a {R n {x) - L n (x))/V2. (15) 



By introducing the Majorana fermions 



Rn = (&„-! + i&n)/V2 



Ho = ^ iv(£ndx£n ~ £nd x £ n ), 



(17) 



n=l 
2N-2 



H' = J2 ~ (18) 



(19) 



In graphene h = m = t±/2 = t/2. The spectrum of the 
Hamiltonian can be derived easily and matches well with 
exact tight-binding results^. 

We emphasize here that the origin of the mass term in 
(fl~9| is due to the fact that in zigzag ribbons an A atom 
from each chain n is connected to a B atom in the next 
chain n + 1 but the B atom in chain n is not connected 
to an A atom in the next chain. In the case of armchair 
ribbons, on the other hand, the connection between A 
and B atoms of adjacent chains is symmetric. However 
this argument does not exclude the presence of the mass 
term in the effective Hamiltonian of armchair ribbons 
because the spinor wavefunction of these ribbons has four 
components and contains two A and two B atoms. 



III. THE SQUARE LATTICE MODEL 



The generalized form of Eq. (|19p with h ^ m, repre- 
sents the square lattice model in Fig. ((T|) which has extra 
hopping terms between B atoms in any given chain and 
A atoms in the next chain. The m (mass) term rep- 
resents the difference between successive rung hopping 
terms, m = tj_(l — 77) /2, and the magnetic term h repre- 
sents the average h = t±(l+rj) /2. Thus, the Hamiltonian 
in Eq. (fi"9"]) describes a continuous family of models that 
include the standard square lattice when 77 = 1 and the ir- 
flux model when rj = — l 17 i 18 i 19 . For the particular case of 
graphene, values of r/ 7^ can represent various hopping 
terms. In particular the NNNN hopping term that has 
been estimated to be of the order of t$ = 0.1 — 0.3e V 20 ' 21 , 
can be incorporated in the renormalization of the mass 
and magnetic terms producing m = (tj_ — 3ts)/2 and 
h — (t± — ts)/2 respectively. Notice also, that this hop- 
ping mechanism is expected to be highly affected by local 
lattice distortions produced either by ripples or bending 
of the material. 

For the generalized case with h 7^ to, following the 
steps outlined in Sec. HI we can solve the tight binding 
Hamiltonian and its wavefunction in a straightforward 



H = 



ip' 
(p 1 



= c 



u A 
u B 



e ia'/2 
± e -i<*'/2 



(20) 



with ip'(k x ,k y ) = t ± (c lk y b + i]c lk « b ) + 2£cos^ and 
ip'(k x ,k y ) = ip'(k x ,—ky). The angle a' is defined by 
(16) ip'(k x , k y ) = ^/(p'ip'c la . The energy of the model is given 
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by E = ±J (tcoa ^ + 2h cos k y b) 2 + 4m 2 sin 2 k y b. As 
before k x is a good quantum number while k y depends on 
k x and is determined by the boundary conditions, thus 
it is a function of W and k x . In the presence of the 
r]t± hopping term the boundary conditions have to be 
generalized so that both components of the spinor wave- 
function are zero at the edges. The wave-function inside 
the ribbon is constructed from the superposition of four 
degenerate states at k y = (±fci,±fc 2 ): 

k y =±k 1 ,±k 2 ^ ' 

6 ^ -<p(-y) ) ° 

After applying the boundary conditions we get 

fk 2 W-a 2 +5\ . (2k 1 y + a 1 +8\ 
rty) = sm J sm ^ j 

(k{W -a x + 8\ . /2k 2 y + a 2 + 5\ 
' Sm { 2 J Sm { 2 J (21) 

with 5 = 0, 7r. Here fci and fc 2 are calculated by the 
degeneracy condition E(k{) = E(k 2 ) and 

Sm { 2 J Sm (, 2 J = 

. fk 2 W-a 2 + 5\ . fk 1 W + a 1 +5\ 
sin I Ism I I (22) 

These two conditions determine k\ and k 2 as functions 
of k x and W . In general fci and k 2 are complex numbers. 

As mentioned above, the limit rj = — 1 describes the 
7T— flux phase that in this language implies the case with 
h = 0. In this regime, the wavenumbers ki and k 2 are 
independent of k x and are given by fci 6 = tt/2 + ri7r6/VK 
and fc 2 & = 7r — fci 6. The energy takes the simple form: 

= ±y^ 2 cos ^ 2 + sin 2 k y b. The valence and con- 
duction bands coincide at four points of the Brilliouin 
zone {k x a = ±TT,k y b = ±7r). At each point there are 
two gapless modes with E = ±2£cos^p w ±tk. Note 
that for odd N one of the solutions of k y lies exactly at 
the gapless corner of Brillouin zone. On the contrary for 
even the zero energy points are always avoided. By in- 
creasing 77, the wave numbers fci and k 2 become complex 
quantities and acquire a dependence on k x . However, for 
odd N , the zero energy mode survives all along the line 
— 1 < i) < 1 with a dispersion that changes from linear at 
rj = -1 to Eq. ^ at T) = 0. As Fig. © shows, the wave 
functions of different modes shows incipient localization 
as 77 increases and become localized when 77 — > 0. In the 
regime with 77 > 0, the magnetic term dominates and 
introduces new zero-energy points for both values of A, 
odd or even. However, for odd N, there is an odd num- 
ber of zero-energy points, one of them being the existing 
mode at k x a = it. In contrast, for even A, there is an 




FIG. 2: The plot of pa— Pb as a function of y across the ribbon 
with W = 46. As approaching the 77 = point (graphene 
ribbon) the wave-function becomes more localized. 

even number of zero-energy points with no zero modes at 
k x a = 7r. Note that the extreme limit 77 — > 1 represents 
the simple square lattice where k y b = nirb/W and energy 
spectrum E = ±(/;cos + tj_ cos k y b). In this case, the 
mode with linear dispersion at k x a = 7r exists only for 
ribbons with even W/b (or odd A). 

In order to visualize the evolution of the band struc- 
tures and wavefunctions of all these different models, it is 
instructive to study the simplest examples of ZGRs with 
A = 2 (even) and N = 3 (odd) number of chains. 

The Hamiltonian of N = 2 ZGR consists of two inde- 
pendent sets of Majorana pairs, odd and even: (£1,^3) 
and (£2, £4)- Consider only one of these sectors: 

+*m(£if 3 - £16) (23) 
By the canonical transformation 

fi = (Xi + X2VV2 £3 = (Xi - X2VV2 (24) 
the above Hamiltonian transforms to 

H = H + im(xiXi ~ X2X2) 

+ih(xiX2 ~ X1X2) (25) 

At 77 = — l(h = 0) the spectrum consists of two de- 
generate massive fermions with E = ±\^m 2 + k 2 . For 
generic 77 the magnetic term h lifts the degeneracy and 
E = ±(v 77i 2 + k 2 ±h). Graphene ribbons are represented 
by the critical value h = m with the particular dispersion 
of edge states given by: E w k 2 /2m. For 77 > the bands 
cut the Fermi energy at two points around k x a = n. 

For ZGRs with A = 3 the Hamiltonian (flU|) consists of 
two independent sets of three Majorana fermions (even 
and odd). Consider one of these sectors: 

H = H + im{^ 3 + &f 5 - 66 - 66) 

-*X66 + 66 - 66 - 66) (26) 
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It is useful to define linear combinations of the Majorana 
fields & and £ 5 as: 7] X = (& + £ 5 )/\/2 and ?] 3 = - 



^5)/v2. Working on the basis 771,773 and £3 it is easy to 
see that the mass term does not couple 773 to the two other 
modes rji and £3 and thus, this mode remains massless. 
By applying a second transformation: 

xi = (6 + m)/V2 xi = (Is + m)/V2 

X2 = (6 - 771VV2 x2 = (&-m)/V2 (27) 

The Hamiltonian transforms into: 

H = H +imV2(xiXi~X2X2)+ih(xiV3+X2V3-XiV3-X2V3) 

(28) 

At 77 = —l(h = 0) the spectrum has two linear modes 
E = ±fc and two copies of massive modes with disper- 
sion E = ±y/k 2 + 2m 2 . The linear modes correspond to 
the k y b = ix modes of the 7r— flux lattice model. An 
increase of 77, introduces the magnetic term that cou- 
ples these modes to the massive modes. As a result, the 
dispersion of the massless modes changes, however they 
remain massless. 

Finally, the energy spectrum of the model for a generic 
value of 77 is given by 

E(E 2 -k 2 ~2m 2 ~2h 2 ) = ±k(E 2 ~k 2 -2m 2 + 2h 2 ). (29) 

When h = m (at 77 = 0) the massless modes dispersion 
becomes E ~ ±fc 3 /4m 2 which describe the edge state of 
the corresponding ZGR ribbon. In the regime 77 > 1 two 
extra zero-energy points are introduced. As 77 — > the 
dispersion of the massless modes becomes linear again. 
These linear modes are in fact the solution of the square 
lattice spectrum near the corners of the Brilliouin zone 
(k x a — ±7r, kyb = ±7r/2). 



IV. EXAMPLES OF SINGLE-PARTICLE 
INTERACTION TERMS 

Another useful aspect of the coupled chains representa- 
tion is that it allows to study the effect of different single- 
particle terms in the Hamiltonian and to obtain the cor- 
responding band-structures in a rather straightforward 
manner. To exemplify these points, in this section we 
derive the Majorana representation of three such terms: 
chemical potential, second neighbor hopping term and 
spin-orbit interactions. We also present the correspond- 
ing expressions for terms that involve spin-orbit inter- 
actions, focusing on the intrinsic and Rashba spin-orbit 
terms. 

It is interesting to notice that the mapping to Majo- 
rana fermions, preserves the quasi-degenerate features of 
the edge bands while captures the main effects of these 
terms in the model for ZGRs. At the same time, the 
total Hamiltonian becomes richer and, as we discuss in 
the next section, some of these terms have direct inter- 
pretations in terms of models of quantum spin chains. 



From Eq. (JT5J) we can obtain the relation between 
the fermion lattice operators c^/b and the Majorana 
fermions representation as: 

c A n (X) = I -1 ) o 



The reverse relations read: 

6n-l +£2n-l = i {- 1 Y{c\ n { X ) ~ C A n (x)) 
6a» + 6ta = -(-l) X (c f An (x) +C An (x)) 
J 





2 






H(6n- 


-1 + £,2n-l, 




+ &»)] 




2 






[(£2>i-l 


~ &n-l) + 


»(6n - 


&»)] 




2 






[fen-l 


— £271-1) — 


*(6n ~ 


Can)] (, 



6n-l ~ &n-l = (-l) :r (ct (x) + C Bn (x)) 



&„ - £a„ = l(-lf(J Bn (x) - C Bn (x)) (31) 

A chemical potential term is given by Ti^(x) = 
pp+{x) = n(pA n {x) + pB n (x)). In the Majorana fermion 
language it reads: 



T~tfj, ~ i2fJ,(^2n-l£,2n + £,2n-l£,2n) 



(32) 



The effect of this term is trivial: it shifts the Fermi energy 
from Ep = to Ep = p, as expected. 

A second-neighbor hopping (NN) term is given by: 

Hnn = <2X)« 3 >( c A i ( x ) c ^(a ; ) + c b,(x)cb j (x)) + h.c. 
The coupling constant £2 has been estimated to take its 
value in the range ti = 0.02/; — 0.2rj eV— . In terms of 
Majorana fermions, this term introduces a velocity renor- 
malization (which can be ignored) due to the inter-chain 
coupling and an intra-chain contribution given by: 



H 



NN 



— l2t2{£,2n-l£,2n + £2n-l£lra) 



(33) 



This expression suggests that the NN hopping term, in a 
first approximation, acts as an effective chemical poten- 
tial, thus breaking particle-hole symmetr y 2 ^! The edge 
state energy is lowered with the corresponding modes ac- 
quiring a finite velocity along the ribbon and becoming 
more stable Fig.©. 

A staggered sublattice potential is represented by a 
Hamiltonian TL v (x) = vp-(x) = v{pA n {x) — pB n {x)). In 
the Majorana fermion representation it takes the form: 



H v {x) « i2l/(£ 2 n-l£2n +£2n-l6n) 



(34) 



For the ZGR in terms of coupled chains, this term opens 
a gap in each single chain (even when decoupled from 
the rest, with t± =0). Thus, every chain making the 
ribbon is gapped or equivalently, each pair of Majorana 
fermions is massive. The total spectrum of the ribbon 
becomes massive regardless of N being odd or even 
Note that this term corresponds to a large momentum 
transfer between two in-equivalent Dirac points. 
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FIG. 3: Energy bands of the zigzag ribbon with N — 3 in the 
presence (solid lines) and absence (dotted lines) of the next 
nearest hoping term t-z = O.lt. 
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FIG. 4: Energy bands of the zigzag ribbon with N = 3 in the 
presence (solid lines) and absence (dotted lines) of the next 
staggered sublattice potential term v — 0.2t. 



Notice that, as a consequence, in iV-odd ribbons, the 
only term which can open a gap at k x a = it, involves 
inter-valley scattering. Although in graphcne the stag- 
gered chemical potential is expected to be zero, recent 
experimental work on suspended graphene samples sug- 
gests the presence of a small energy difference between 
the two sub-lattices^. Moreover in nanoribbons made 
of SiC or BN there is an intrinsic staggered chemical 
potential due to the presence of two different atoms in 
sublattices A and BM£Z£& 

The procedure can also be applied to other single- 
particle terms in the Hamiltonian. For example, we 
can also study the expression of spin-orbit interaction 
terms using the Majorana representation. The intrinsic 

spin-orbit (I-SO) is Hjso = H' s Yl<ij> v ij c t c j m which 
s = ± refers to the z-component of the electron spin 



operator. Vij = 2(dik x dkj)/V3 is the path-dependent 
number which takes the values Vij = (±1) depending on 
the shortest path for a hopping process that goes from 
site i to the next nearest cite j through an intermediate 
cite fc 29 i 30 ' 33 . For modes with spin up the spin-orbit is 
written as 

Hi SO ~ i'(£nfn+2 + (35) 

For spin down the equivalent term corresponds to the 
transformation t' — > — t'. It has been shown that in the 
presence of the I-SO interaction, the bulk states have 
a gap of the order of t'. The situation is different in 
ribbons. For ZGRs with N odd, the edge states get ex- 
changed as the conduction and valence bands cross each 
other at the midgap point k x a = ir. The ribbon thus re- 
mains metallic. In this case the I-SO interaction removes 
the degeneracy of the edge states near k x a = ir and ren- 
ders a pair of linear dispersion bands. Moreover, the edge 
states become spin-filtered with opposite spin currents at 
opposite edges of the ribbon. This is the topological in- 
sulating or quantum spin Hall phase of graphene 2 ^. For 
ZGRs with even N the ISO interaction opens a gap be- 
tween the edge states at k x a = ir—. The gap is very 
small and scales as A ~ (t'/t) N . Graphene ribbons in 
this phase are normal insulators Fig. ([S"|). 

For graphene samples on a substrate or suspended but 
subjected to a perpendicular electric field, the Rashba 
spin-orbit (RSO) interaction becomes quite relevant. 
The expression for the RSO interaction is Hrso = 
i'\Y. where £ z . Here dij is a vec- 

tor which connect the nearest neighbor cites i to j. The 
RSO has two contributions: the intra-chain part gives 

HrSOI = — *4v3A(£ 2 nT£2ni — &n^2nl 
+£2n-lT£2n-lJ. — ^2n-lT^2n-lJ.) (36) 

The inter-chain part is 

HrS02 = i2A(^2n-lT^2n+21 - ^2n-lT^2n+2J. 
+&n-ll£,2n+11 — ^2n-lJ.^2n+lT 

— C2nT?2n+l| + ^2nT^2n+l| 

— C2ni?2n+1T + ^2ni^2n+lT 

— £,2n- ll£2n+2 1 + ^2n-l|'?2n+2J. 

— £2n-ll£2n+2T + £,2n-l 1&TI+21 
+C2nT^2n+U — C2nT^2n+ll 

+^2ni^2n+lT — ?2nl^2n+lT) (37) 

Notice that the RSO does not open a gap in the bulk^i. 
The same holds for ZGRs: there is no gap in the ribbon' 
spectrum and the dispersion remains power-law like near 
the zero-energy point. However, in the presence of RSO 
interactions, the zero-energy point located at k x a = tt 
for RSO zero, shifts in opposite directions for opposite 
spins 3 ^ Fig. ((6]). For wider ribbons with quasi- flat (and 
originally spin-degenerate) bands, there is a more pro- 
nounced lift of the spin-degeneracy. Another interesting 
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FIG. 5: Energy bands of the zigzag ribbon with N = 2 (top) 
andiV = 3 (bottom) in the presence of the intrinsic spin-orbit 
t = Q.lt. For even N there is a gap at k 



point is that the RSO interaction does not remove the 
even-odd width dependence of the edge states bands for 
finite widht ZGRs. As a consequence even in the pres- 
ence of RSO a ribbon with odd number of chains re- 
mains metallic while one with even number of chains is 
insulating. Notice also that these results have important 
consequences for the analysis of the topological insulator 
phases that can be observed in graphene ribbons as dis- 
cussed in Refi^. This can be readily seen by considering 
an even N ribbon: in this case, there is an even number 
of Kramer pairs at each edge. In contrast, for odd N 
ribbons, the number of Kramer pairs is always odd, even 
in the presence of RSO interactions. 



V. GRAPHENE RIBBONS AND QUANTUM 
SPIN CHAINS MODELS 

The Hamiltonian (fl^ can also describe some specific 
quantum spin chains models. Before presenting the con- 
nection between the two approaches, it is instructive to 



k/7T 

FIG. 6: Energy bands of the zigzag ribbon with N — 3 in the 
presence (solid lines) and absence (dotted lines) of the Rashba 
spin orbit interaction A = O.lt. 

briefly review some results of the models of: two coupled 
spin-1/2 chains and a spin-1 chain^ 
a) Two coupled spin-1/2 chains. 

The symmetry of two spin-1/2 Hciscnbcrg antiferro- 
magnetic chains, in the absence of the inter-chain cou- 
pling, is SU(2) x SU(2). In the continuum limit the 
spin operator in each chain, can be written in terms 
of current operators and a staggered part as Si(x) = 
Ji(x) + Ji(x) + (— \) x ' a rii{x) where i = 1,2 labels chain 
i. The smooth part of the spin operator corresponds 
to the sum of right and left mover currents (J + J) 
of the SU(2) k=1 WZNW model, and represents the to- 
tal magnetization of the chain. The staggered magneti- 
zation is represented by n and generalized models con- 
tain a dimerization operator e = ( — l) x ' a S(x).S(x + a) 
that represents the time component of the staggered 
magnetization^. 

The addition of an interchain coupling reduces the 
symmetry to the generic SU{2) x Z2. In the standard 
model for two spin-1/2 coupled chains (ladder)^, the in- 
terchain coupling is isotropic: H' = J1S1.S2. The rele- 
vant part of this term stems from the staggered magneti- 
zation contribution and is given by H' as J " \ri\.rii- The 
consequence of including such a term is the opening of a 
gap in the spectrum. In terms of Majorana fermions, the 
total Hamiltonian can be written in terms of four fcrmion 
fields £1 , £2 , £3 and £4 . Three of these fields £ = (£1 , £2 j £3) 
have equal masses m ~ Jj_ and the remaining one (£4) 
has a mass —3m. The classification of Majorana fermion 
states into triplet and singlet states reflects the underly- 
ing global SU{2) symmetry of the ladder. 

Majorana fermions £ describe the triplet excitations of 
the two-coupled 5=1/2 quantum chain^i. When the 
model is written in terms of right and left movers it can 
be shown that the total currents I = J\+ J2 for each type 
of mover satisfy the Kac-Moody algebra of the SUk=2(2) 
vector currents. Since the central charge of the model is 
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C = 3/2, the currents can be represented using a triplet 
of Majorana fermions: I a = -^e abc £h£ c where a, b, c refers 
to 1, 2, 3. The model has two primary fields namely tj^ 1 / 2 ) 
and with the second field given by $^ = £ a £b- The 
other Majorana field £4 describes the singlet excitations 
of the model. This field £4 contributes to the relative (or 
the so called wrong^i) currents K = J\ — J 2 which are 
represented by K a = i^a^i- It is convenient at this point 
to introduce the fields Q a = £4^ and Q a = ^£, a . 

Notice that the ferro-antiferromagnetic transformation 
J± — ► — J±_ (which is reflected in to — ► — m) does not 
change the thermodynamics but it interchanges singlet 
and triplet roles. 

b) Spin S — 1 chain 

The connection to a S = 1 Heisenberg spin chain is 
given by considering the biquadratic S = 1 chain as de- 
fined by2&29 

H = J2[S(x).S(x + a) + b(S(x).S(x + l)) 2 ] (38) 

X 

At b = — 1 the model has gapless modes and it is equiv- 
alent to a WZW SU(2)k=2 model. The Hamiltonian at 
this point can be described by Majorana fermions as 

3 

Hq = ^2 iv (£nd x £,„ - £rA£n) (39) 
ri=l 

For b < — 1 there are two degenerate massive dimerized 
phases. For — 1 < b < 1 the spin chain is in the Haldanc 
phase with a gapped singlet ground stat o 40 i 41 . 

The link between the spin operator and the WZW 
model is given by 

S(x) ~ I(x) + I{x) + {-l) x n(x) (40) 

Here n is the staggered part of the spin field. In analogy 
with the coupled spin-1/2 chains, the smooth (average) 
magnetization is given by M = 1 + 1. 

For < b < — 1 the model has a gap that increases 
monotonically as a function of to = —1 — 6. Close 
to b = — 1 the model can be described effectively by 
adding relevant mass terms. The mass term is positive in 
the dimerized phase and negative in the Haldane phase. 
These mass terms are expressed in terms of the relevant 
fields = £ a £b with a 3 x 3 matrix representation. For 
isotropic models the mass term is written by mTr($ 1 ). 

c) ZGRs and quantum spin chains 

To show the connections between these models and 
graphene ribbons we start with the even N = 2 ZGR. 

The Hamiltonian written in terms of right and left 
movers can be bosonized in a straightforward manner 
by introducing the bosonic fields <ft+ = exp(±iV^l7r</> + ) = 

(£i + *&)/\/2 and 4>+ = exp(±»V^+) = + i&)/V2 
for one set of movers and the corresponding 0_ , 0_ for 
the other. The total Hamiltonian density reads 

H = Ho + i— V sin V4^9„ + -A= V vd x Q v (41) 



where 6„ = 0„ — <f> v is the dual field of <f>„ = <p v + 
4> v . Equation (|4T|) describes two decoupled sine-Gordon 
Hamiltonianas with a magnetic field (or chemical poten- 
tial) applied to the dual sector. 

The model for the N = 2 ribbon can be mapped to 
the two-leg ladder Hamiltonian discussed above with an 
anisotropic interchain coupling and an in-plane applied 
magnetic field: H = H + H' lh + H' lm + H' 2m . The 
magnetic field term is expressed in terms of the current 
operators as H' lh = h(K l J - K y + I v - I v ) that are 
written in terms of Majorana fermions as shown above. 
The operator content of these currents can be given in 
terms of vector spin-chirality operators as: J — I ~ 
Si{x) x Si{x + a) + S 2 (x) x S 2 (x + a) and K - K ~ 
S\{x) x S\(x + a) — S 2 {x) x S 2 (x + a). The mass terms 
are produced by the following two interchain couplings. 
The first term is Ti' lm = im(Q 2 — Q 2 ). Expressed as spin 
operators this is Q 2 — Q 2 ~ (Si x S 2 ) v ~ (ni x n 2 ) v . 
The second term is 7i.' 2m = im(^^ + Q^i) which can be 

generated from (3? 13 + $31) ~ i n i n 2 + n 2 nf). 

The analysis can be extended to the next even ZGR 
ribbon, namely the N — 4 ribbon which can be described 
by two sets of four Majorana fermions. The Hamiltonian 
of each sector, separately, can be represented by a two- 
chain spin-1/2 ladder. To get the equivalent spin chain 
model, we just need to add two more terms to the equiva- 
lent spin model for the N — 2 ZGR discussed above: The 
first term is H' 3m = im(Qx - Qi) - t±(Si x S 2 ) x . The 
second term is the magnetic field H 2h = h(K x — K x ). 

The mapping between odd ZGR ribbons and spin- 
chain models is simpler for the N = 3 case. The equiv- 
alent spin-chain describes the anisotropic version of the 
biquadratic spin S — 1 chain2& with an in-plane magnetic 
field. The Hamiltonian of one of the sectors is given in 
(|2"5)l . This model describes aS=l model with masses 
mi = to, 7Ti 2 = —to and to 3 = and two magnetic fields 
hi = h 2 . In the language of ZGR ribbons, we notice that 
the value of the mass terms mi and m 2 can be changed 
when including the ISO interaction (|35[) . In this case 
then, the mass terms are given by mi = t±/2 + st' and 
to 2 = —tj_/2 + st' where s = ± stands for real spin-up 
spin-down electron. 

For the TV = 5 ZGR each sector has five Majorana 
fermions. In principle the Hamiltonian can be described 
by the S = 2 representation of the SU k=w {2) WZW 
model. The corresponding primary field has the confor- 
mal dimension (1/2,1/2) which can be written in term 
of Majorana fermions bilincars. Furthermore, the central 
charge of the theory is C — 5/2 which is another indica- 
tion that five Majorana fermions are needed to describe 
the model. 



VI. CONCLUSIONS 

Zigzag graphene ribbons show remarkable physical 
characteristics that are predicted to have important ef- 
fects in their transport properties. In this work we have 
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provided a detailed analytic treatment of a tight-binding 
Hamiltonian with hard- wall boundary conditions that ex- 
plains the surprising width-dependent properties in terms 
of the existence of a zero-energy mode that corresponds 
to localized states along the edges of a finite width rib- 
bon. As a consequence, ZGRs with even widths (odd 
number of chains) are metallic while odd width ribbons 
(even number of chains) are insulating. 

The unusual even-odd dependence led us to consider 
models of ZGR in terms of coupled quantum chains. This 
particular continuous limit preserves the main features of 
band structures and spinor wave-functions while keeping 
the width-dependence. 

In the Majorana fermion representation, the model of 
ZGR in terms of coupled chains can be easily extended 
to include various single-particle interactions and we have 
analyzed in details the effects of adding a chemical poten- 
tial term, second and third nearest neighbor hoppings, 
a staggered chemical potential, intrinsic spin-orbit and 
Rashba spin-orbit interactions. 

Besides the straight-forward treatment of these terms, 
the mapping reveals that ZGRs can be viewed as a mem- 
ber of a continuous family of models in square lattices 
that range from the standard square lattice to it— flux 
models. These models share the same even-odd width 
dependence first obtained in ZGR models. We have pre- 
sented a full solution for these models and emphasized 



their similarities with ZGRs. 

As a final application, we have shown that it is pos- 
sible to treat ZGRs as coupled quantum spin chains. A 
careful analysis of models in terms of Majorana fermions 
shows that ZGR with two chains inside are fully equiva- 
lent to two spin- 1/2 coupled chains with an anisotropic 
interchain coupling and an external magnetic field. We 
also analyzed another example of an even-width ribbon 
(N = 4). We applied the mapping to odd- number mod- 
els and showed that an N = 3 ribbon, corresponds to the 
anisotropic biquadratic spin 5 = 1 chain with an in-plane 
magnetic field. A preliminary analysis of an TV = 5-chain 
ribbon suggests a description in terms of a set of 5 Majo- 
rana fields for each Hamiltonian sector. We would like to 
remark that the mapping proposed here allows to obtain 
solutions for all these new quantum spin chains models 
by direct comparison with the known results obtained 
with ZGR's. 
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